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Abstract 

It is shown, with the example of the experimentally known Adler function, that there 
is no matching in the intermediate region between the two asymptotic regimes described 
by perturbative QCD (for the very short-distances) and by chiral perturbation theory 
(for the very long-distances). We then propose to consider an approximation of large-iVc 
QCD which consists in restricting the hadronic spectrum in the channels with quan- 
tum numbers 0~, 0+ and 1+ to the lightest state and to treat the rest of the narrow 
states as a perturbative QCD continuum; the onset of this continuum being fixed by 
consistency constraints from the operator product expansion. We show how to construct 
the low-energy effective Lagrangian which describes this approximation. The number of 
free parameters in the resulting effective Lagrangian can be reduced, in the chiral limit 
where the light quark masses are set to zero, to just one mass scale and one dimensionless 
constant to all orders in chiral perturbation theory. A comparison of the correspond- 
ing predictions, to 0{p'^) in the chiral expansion, with the phenomenologically known 
couplings is also made. 



1 Introduction 



We shall be concerned with Green's functions of colour singlet operators which are bilinear in 
the light quark fields. There are at present two asymptotic regimes where QCD is a predictive 
theory for these Green's functions: the short-distance regime where the external momenta 
are deep euclidean in regions of non-exceptional momenta |j^], and the very long-distance 
regime where the external momenta are very soft. The behaviour of the Green's functions in 
the short-distance regime can be successfully predicted via perturbation theory in powers of 
the QCD running coupling constant (pQCD). In the chiral limit where the light quark masses 
are set to zero, the behaviour in the short-distance regime is then governed by only one mass 
scale, e.g. the parameter for three flavours. 

The way that QCD predicts the behaviour of the Green's functions in the long-distance 
regime is not so straightforward. At very low energies, the sole input of the property of 
spontaneous breakdown of the chiral SUl{3) x SU{3)r invariance of the QCD Lagrangian with 
massless u, d, s quarks, allows one to set up a well defined theoretical framework |2|-|6| known 
as Chiral Perturbation Theory (xPT). This is formulated in terms of an effective Lagrangian 
that describes, at low energies, the strong interactions of the Goldstone modes associated 
with spontaneous chiral symmetry breaking. The behaviour of the Green's functions in the 
very long-distance regime is then predicted by this effective Lagrangian as an expansion in 
powers of momenta. This, however, comes at the expense of introducing an increasing number 
of local operators of higher dimension with unknown couplings which build the effective 
Lagrangian at higher orders in the chiral expansion. [Some of these local operators act as 
counterterms in the loop expansion generated by the lower order terms.] At 0{p^) in the chiral 
expansion the number of independent local operators is already of the order of a hundred Q. 
Although in principle their couplings are fixed from QCD, their calculation demands a genuine 
knowledge of non-perturbative QCD properties. In certain circumstances one can relate these 
low-energy constants to a sufficient amount of independent experimental observables which 
determine them and many successful predictions have been made this way |^. [This is the 
orthodox methodology of the science.] However, as higher orders in the chiral expansion 
are required, it often happens that there is simply no model independent way to fix all the 
couplings. In practice the onset of this situation is at 0{p^) in the chiral Lagrangian for the 
strong interactions and already at ©(p^) in the sector of the non-leptonic weak interactions. 
If one wishes to make progress here, one clearly needs to find a good theoretical approximation 
substitute for full fledged QCD. 

An interesting question one may ask is: what is the matching of the two asymptotic 
regimes which we have just described? This is precisely the question which we wish to discuss 
in this paper. For simplicity we shall illustrate most of this discussion with a specific Green's 
function, the so called Adler function (see section 2 below for definitions), but many of the 
features presented here are rather generic. We claim that the direct matching is in fact 



extremely poor; even when non-perturbative short-distance power corrections a la SVZ [11| 
are retained. Once more we find that if one wishes to make progress in the quality of the 
matching between the pQCD short-distance behaviour and the behaviour which is only 
operational at very long-distances, one clearly needs a good theoretical approximation as a 
substitute for full fledged QCD. 

As a first step in this direction we shall consider QCD in the limit of a large number of 



^See e.g. the review articles in refs. 
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colours |T3|-[l7|. This limit, which following 't Hooft's notation we shall term QCD(cxd) for 
short, offers a convenient starting point where one can hope to discuss fruitfully the afore- 
mentioned matters since it is still rich enough to share with QCD its most relevant non- 
perturbative features. Among them, a very important one, is the fact that QCD(oo) has the 
correct pattern of spontaneous chiral symmetry breaking |[l^ , p^ , provided that the property 
of confinement still persists in going from QCD to QCD(oo). The advantage of the large— A'c 
limit is the emergence of narrow meson states in the Green's functions. The disadvantage is 
that the number of narrow states is infinite and they appear at all energy scales, even when 
the euclidean momenta are large enough that a perturbative treatment should be adequate. 
Obviously, a perturbative description, when possible, is more predictive than an infinite tower 
of narrow states with masses and decay constants which are a priori unknown. 

We propose to consider an approximation to QCD(oo) which consists in restricting the 
hadronic spectrum in the channels with quantum numbers 0~, 1^, 0^ and 1^ to the 
lowest energy state and in treating the rest of the narrow states as a QCD(cxd) perturbative 
continuum; the onset of the continuum being fixed by consistency constraints from the op- 
erator product expansion (OPE) as we shall discuss. The reason why we limit ourselves to 
this simple hadronic spectrum is that, in the large- A'^c limit, this is the minimal structure 
compatible with spontaneous chiral symmetry breaking The traditional successes of 

"vector meson dominance" in describing low energy hadron phenomenology guarantees that 
this approximation, which we shall call "lowest meson dominance (LMD) approximation to 
large- Ac QCD", should be rather reasonable 0. The immediate question which then arises 
is the following: is there an effective Lagrangian formulation which explicitly exhibits the 
properties of QCD(oo) in the LMD approximation? We shall show that there is indeed such 
an effective Lagrangian. The effective Lagrangian in question can be obtained starting from 
a modified version of an extended Nambu-Jona-Lasinio type Lagrangian (ENJL) which has 
been very much discussed in the literature The modifications, however, are highly non 
trivial because they require the presence of an infinity of local operators which, however, have 
known couplings. The couplings of the new operators are known because they are fixed by 
the requirement that they remove, order by order in the chiral expansion, the effects of non- 
confining QQ discontinuities produced by the initial ENJL ansatz. The resulting Lagrangian 
when the quark fields are integrated out is a non-linear sigma type Lagrangian with degrees 
of freedom corresponding to the Goldstone fields and also to Vector, Axial- Vector, and Scalar 
fields. The matching of the two-point functions of this effective Lagrangian to their QCD 
short-distance behaviour can be systematically implemented. However, for Green's functions 
beyond two-point functions, the removal of the non-confining QQ discontinuities produced 
by the initial ENJL ansatz is not enough to guarantee in general the correct matching to 
the leading QCD short-distance behaviour and further local operators have to be included. 
In the sector of Vector and Axial-Vector couplings, and after the matching with the QCD 
short-distance behaviour from the OPE for three-point functions is implemented |^6|, the 
resulting low-energy Lagrangian to 0{p^) in the chiral expansion coincides with the class of 
Lagrangians discussed in ref. |27]. The advantage, however, with respect to a purely phe- 
nomenological description in terms of chiral effective Lagrangians which include resonances 

^If necessary, one could also add further states, (see e.g. ref. for a recent discussion of a possible 
description of the 7r(1300).) However, as we shall see, the phenomenological success of the LMD approximation 
indicates that, to a good approximation, these higher states play a minor role in the dynamics of the low-energy 
effective Lagrangian. 

''See e.g. refs. |2ll-M] and references therein 
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as discussed e.g. in refs. [p8|-|31| and references therein, is that in the LMD approximation 
to QCD(oo) , the number of free parameters can be reduced to just one mass scale and one 
dimensionless parameter to all orders in xPT. 

2 The Adler Function 

We are concerned with two-point functions 

n^.(g)a6 = i J d^xe^'i-^ < o|r{y;(x)K'(o)}|o > (2.1) 

of vector quark currents 

V^{x)=q{xh^^qix), (2.2) 

where A" are Gell-Mann matrices (trA"A'' = 26ab) acting on the flavour triplet of u, d, s light 
quarks. In the chiral limit where the light quark masses are set to zero, these two-point 
functions depend only on one invariant function (Q^ = —q^ > for spacelike) 

^fiu{q)ab = {qtJLqv - g^iuq^)^{Q^)^ah ■ (2.3) 

In QCD, the invariant function Il{Q^) obeys a once-subtracted dispersion relation and it is 
conventional to choose the subtraction at = 0: 

n(Q2) = n(o) - q2 / ^^^-imn(t). (2.4) 

The Adler function is defined as the derivative of the invariant function n(Q^) i.e., 

-A(Q^) = -Q^^^, (2.5) 

and therefore, its relation to the spectral function ;ilmn(t) does not depend on the choice of 
the subtraction 

/•oo q2 1 



2.1 Short— Distance Behaviour 

The Adler function obeys a homogeneous renormalization group equation, and it has been 
computed in pQCD (in the MS renormalization scheme) up to ©(a^), with the result |34,35| 
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and Uf the number of hght quark flavors. The behaviour of the Adler function as a function 
of predicted by pQCD is shown in Fig. 1. The plot in Fig. 1 corresponds to the choice 



Q in eq. (2^), with as{Q ) the solution to the implicit equation 
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which results from solving the running coupling constant renormalization group equation 
at the two-loop level. Figure 1 clearly shows the asymptotic freedom behaviour: A{Q^ — > 
as well as the fact that at < IGeV^ we are already entering a region where 
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non~perturbative effects are likely to be important. 
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Fig. 1 The Adler function from, perturbative QCD. The dashed lines correspond 



to A- 



MS 



(372±76)MeV ||,|7|. 

As first pointed out by Shifman, Vainshtein and Zakharov |11|, non-perturbative 
power corrections appear naturally when two-point functions like eq. (2.1) are considered in 
the physical vacuum within the framework of the operator product expansion (OPE). The 
leading ^ and next -to-leading ^ non-perturbative corrections were calculated in ref. | pT[ . 
They depend on the size of the gluon condensate, which gives the leading correction to the 
perturbative result in eq. (|2.7|): 
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and on the size of four quark condensates, which give the next -to-leading power correction 
to the perturbative result in eq. ( |2.7| ): 



3 ((ti7/.75A"n - d-ff^j^X^'d)' 



q=u,d,s 



vros 



(2.12) 



To leading order in the l/A^c~expansion four quark condensates factorize into products of the 
lowest dimension quark condensate with the result 



2\Nc^oo 



28 vras 

T 



(2.13) 



The phenomenological determinations of the gluon condensate and the quark condensate have 
unfortunately rather big errors. A generous range for the gluon condensate which covers most 
of the determinations is 

(asG^G" ) = (0.08 ± 0.04) GeVl (2.14) 



The quark condensate is scale dependent. The values obtained from a variety of sum rules 

(V^V) (1 GeV^) = [(-235 ± 27) MeV]^ . (2.15) 
Using these values, the resulting prediction for the Adler function is shown in Fig. 2. 
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Fig. 2 The Adler function including (a) the corrections, (h) the and corrections. 
The dashed line is the VMD prediction corresponding to eq. \2.1f\ ) below, valid at low Q. 



See e.g. the recent compilation of sum rule estimates in ref. Our error in eq. (2.15) is however larger 
than the one quoted in this reference. 
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The curve (a) in the figure is the one which results from adding only the ^-corrections to 
the pQCD result shown in Fig. 1; the curve (b) includes both the ^ and ^-corrections. 
The curves in Fig. 2 clearly show that with the values for the lower dimension condensates 
given in eqs. ( 2.14| ) and (2.15), the contribution from the leading and next-to-leading non- 
perturbative power corrections -at values where one can trust the pQCD contributions- 
turn out to be already quite small. This is due to the empirical fact that the scale of per- 
turbation theory is rather large (we are using [^,|3^ ^ms — (372 ± 76) MeV in our 
plots). For values of as large as that, one is forced to choose > 1 GeV^ to trust the 
perturbative expansion in powers of as{Q^), and for > 1 GeV^ the power corrections are 
already quite small. 



2.2 Long— Distance Behaviour 

The behaviour of the invariant function Il{Q^) in eq. ( p. 3D at = is governed by a 
combination of two 0{p'^) coupling constants of the xPT effective Lagrangian: 

n(Q2) = _4(2Fi + Lio) + 0(Q2). (2.16) 
The constants Hi and Lio are the coupling constants of the terms 

C^six) = ■■■ + Li^iTU{x)^F^''{x)U{x)FL^,u{x) 

+Hitr {F>^' {x)Fr^,{x) + Fr(x)Fi^.(x)) + • • • (2.17) 

in the effective chiral Lagrangian, where Fl and Fr are the field strength tensors associated 
with external — and + matrix field sources, and U{x) the unitary matrix which 
collects the Goldstone fields {U VrUVI under chiral transformations). The Adler function 
is not sensitive to n(0), and hence to 2Hi + Liq. The slope of the Adler function at the origin 
corresponds therefore to combinations of couplings in xPT which are already of 0{p^) and 
therefore theoretically unknown a priori. There are also contributions to the slope of the Adler 
function from chiral loops involving lower order couplings. These chiral loop contributions 



are calculable [38|, but they are non-leading in the l/A'^c~expansion and therefore we shall 
ignore them here. 

One can make an "educated" guess of the value of the slope of the Adler function at the 
origin by invoking VMD arguments. In this case it amounts to the assumption that, at low 
energies, the vector spectral function is dominated by the narrow width pole corresponding 
to the /9-meson 



1 
vr 

This leads to the "prediction" 



-lmn(t) ~ 2f^M^p6{t - Ml) . (2.18) 



^(Q'^0)vMD = ^Q' + O(Q^), (2.19) 



and it corresponds to the dashed line starting at Q'^ = which we have plotted in Fig. 2. 

We conclude that the only rigorous knowledge we have from QCD about the long-distance 
behaviour of the Adler function is that it vanishes at = 0, with a slope which is governed 
by 0{p^) terms in and therefore its determination is, at present, model dependent. We 
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wish to emphasize, however, that the property that it vanishes at = is highly non-trivial. 
There is nothing in the pQCD behaviour in eq. ( |2.7D which points to that. Although one can 
find in the literature arguments invoking that Q!s(Q^) might "freeze" to a constant value 
when 0, one should be aware that in the case of the Adler function a naive "freezing 

argument" contradicts |^ the long-distance QCD behaviour encoded in xPT. 



2.3 Empirical Determination of the Adler Function 



The isovector hadronic spectral function has been determined both from e^e^ -annihilation 
measurements and from hadronic r-decays ^ We show in Fig. 3 the compilation from the 
T-data as given in ref. |40|, except for the overall normalization. 




0.5 1 1.5 2 2.5 3 

t [GeV^] 

Fig. 3 The isovector hadronic spectral function from hadronic t - decays. The dashed line 



is the asymptotic freedom prediction 



Air-' 



The normalization in our Fig. 3 is in accordance with the definition of the spectral function 
given in eqs. (2.1) to ( |2.4| ) above. When extracting the Adler function from experimental 
data, one is confronted with the question of matching to the asymptotic QCD continuum 
beyond the region accessible to experiment. The concept of global duality, as explained in 
refs. [^,^, provides the way to do this matching. In practice, it consists in introducing an 

^See the appendix in ref. for further discussion on this point. 

^ For a recent review, where references to the original experiments can be found, see ref. |4C|]. 
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ansatz for the complete spectral function 



^ Imn(t) = ilmn(t)exp. 9{so - t) + hmU{t)^QCj) e{t - sq) , 



vr 



vr 



(2.20) 



where ^Imn(f)exp. denotes the hadronic experimental spectral function and ^lmIl(t)-pQQY) 
the spectral function predicted by pQCD. The onset of the continuum sq is then fixed by 
solving the equation in sq: 



/•so \ rso \ 

dt -Imn(t)exp. = dt -Imn(t)pQci3 



(2.21) 



This equation guarantees that the matching between long distances and short distances in the 
Adler function is consistent with the OPE. Indeed, it can be easily shown that the equation 
in ( |2.21[) follows from the absence of operators of dimension two Q in the OPE of the product 
of current operators which defines the Adler function. 

Using the experimental input corresponding to the r-decay data shown in Fig. 3, and the 
expression R 



so 



1 



dt -Imn(t)pQCD 



Nr 4 



+ 
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as(so)^ ^ 



as(so) 



vr 



vr 



+ 



as(so) 



vr 



for the pQCD r.h.s. in eq. ( p. 21 ) results in the solution 

(1.60 ±0.17) GeV^ 



So 



(2.22) 



(2.23) 



where the error comes from the error on A-p^ only. The resulting Adler function is shown in 
Fig. 4 below. 



^ In the presence of quark masses there is in fact a correction term to eq. (2.21) of C'(^) which we are 
neglecting here. 

* For the purposes we are interested in here, the evaluation of this integral as a contour integral in the 
complex t-plane, as proposed e.g. in ref. does not change significantly the result. 
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Fig. 4 The Adler function from experimental data (full line) and from the VMD param- 



eterization (dashed line) in eq.(2.24)- 



It is interesting to compare the "full experimental" result shown in Fig. 4 with the one obtained 
from the simplest hadronic parameterization one can think of, which is the one corresponding 
to a narrow width vector state -the p{770)- plus the corresponding pQCD continuum i.e., 



1 



TT 



1_ 



Imn(t) = 2f'^M'p6{t - M'^) + -Imn(t)pQCD 0{t - s 



(2.24) 



The onset of the continuum for this parameterization is at Sp = (1.53 it 0.06) GeV^, and 
the corresponding Adler function is also shown in Fig. 4 -the dashed curve- . Although 
point by point the two spectral functions in eqs. (2.20) and (2.24) are very different below 
their continuum onset, the associated Adler functions result in rather similar functions of 
the euclidean variable . It is in this sense that these two spectral functions are globally 
dual MM. 



We wish to summarize the basic features which we have learnt in this section. 

• The direct matching between and pQCD is very poor, regardless of whether or 
not one includes non-perturbative power corrections a la SVZ [0]. Including these 
corrections does not help much because of the empirical fact that A^^ is large. At the 

values above which we can trust pQCD, the non-perturbative power corrections are 
already quite small. 

• As shown in Fig. 2 there is no matching between the naive extrapolation of the linear 

behaviour predicted by the first non-trivial order of xPT, Oip^) in this case, and 
the pQCD prediction with inclusion of the leading and next-to-leading ^ non- 
perturbative power corrections. 
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• The simplest hadronic parameterization which consists of a narrow-width vector state 
plus the corresponding pQCD continuum does remarkably well in reproducing the ex- 
perimental shape of the Adler function. This is very encouraging in view of the LMD 
approximation to QCD(oo) which we have suggested in the Introduction, namely to 
consider that the large-A'^c spectrum consists of a prominent low-energy narrow state 
with the rest of the narrow states lumped together as a QCD(oo) perturbative contin- 
uum. It is this LMD approximation to QCD(oo) which we are now going to further 
explore in the following sections. 



3 Effective Lagrangians and Matching Conditions 



Internal quark loops in Green's functions of colour-singlet quark bilinear operators are sup- 
pressed in the large- A'c limit. Quark degrees of freedom appear only in the external loop. At 
low-energies, the external momentum flowing into the quark loop is small. If one "assumes" 
that the dynamics of the gluons has a characteristic non-perturbative scale > 1 GeV ^ , it 
is natural to think of this quark loop as a low-energy insertion of effective higher-dimensional 
composite quark operators. These operators will be local at low-energy scales below A-^. Fur- 
thermore, since the physics encoded in these effective operators originates from integrating 
out degrees of freedom from A^ to infinity, the QCD U{3)l x U{3)r flavor symmetry must 
also be explicitly realized in the effective composite quark operators. Naively, one expects A^ 



and the sq of subsec. ^ which fixes phenomenologically the onset of the pQCD continuum, 
to be rather similar. 

In perturbation theory one can explicitly see the emergence of these higher dimensional 



composite operators at very large orders in the QCD coupling constant 1 39, 44, 45]. Although 
this is not a proof of the existence of these operators, it is nevertheless reassuring to know that 
there are physical systems which are well described by an analogous mechanism. Indeed, a su- 
perconductor is an example of a physical system where the appearance of higher-dimensional 
fermionic operators in the non-perturbative dynamics is also hinted at by the very large or- 
ders of perturbation theory, (see e.g. ref. [^].) It is this kind of reasoning that leads one to 
consider as a plausible starting ansatz to formulate the LMD approximation to QCD(oo) , an 
effective Lagrangian of the form 



+ 



(x) i-i^D^+i{s -i-i^p) q{x) 
Stt^Gs 



A a, b=navour 



A a, b=navour 



E lirLl'-qi) {<fLl,Qi) + (L - R) 



+ • 



(3.1) 



where Dn 



du - il 



1-75 



'/i — '-'^ '■'-/^ 2 ^'"m 2 > — "^M "A" ' M — 'JpL ~r uifj,, 

matrix source fields and Gs and Gy are dimensionless couplings of 0{N^). Summation over 
colour indices within brackets is understood. This effective Lagrangian is a Nambu-Jona- 
Lasinio type Lagrangian [21| in its extended version 25 1 . For the sake of brevity we 
shall refer henceforth to the set of operators explicitly shown in eq. (|3.1|) as ENJL. The 



H75 



+ ttfj,, s and p are external 



Glueballs, for instance, are expected to have masses > 1 GeV. 
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presence of the two types of four-quark interactions is the minimum required to have non- 
trivial dynamics in channels with the quantum numbers 0^ , 1^ , 0^ , and 1^ . Taken as it 
stands this Lagrangian, at leading order in l/N^. and provided that the coupling Gs becomes 
critical {Gs > 1), yields the wanted pattern of spontaneous chiral symmetry breaking and 
reproduces surprisingly well many of the low energy properties of hadron phenomenology 
p5[ with only three independent parameters: Gs,Gy and A^^. 

The ENJL Lagrangian, however, has basically two very important drawbacks. On the 
one hand, although at = it produces the correct nonet of Goldstone poles in the S- 
matrix that are color singlets and have the right flavor quantum numbers; at higher energies 
it has unconfined quarks, whose most evident effect is to create unphysical imaginary parts 
in Green's functions. Furthermore, the unconfined quark-antiquark pairs produce S-matrix 
elements which do not obey the mesonic large- A'c counting rules. On the other hand, there is 
no explanation as yet on how one is supposed to do the matching of this effective Lagrangian 
with short-distance QCD. As we shall see, the two drawbacks are very much related to each 
other and to the ellipsis in eq. (|3.l|). 



3.1 The Adler Function of the ENJL— Lagrangian. 

The explicit form of the Adler function predicted by the ENJL Lagrangian has been calculated 
in ref. [24|. It reads 



•^enjl(Q ) 



(3.2) 



where n(Q^) is given by a loop of constituent free quarks. Since this function is ultraviolet 
divergent it requires the introduction of a cut-off. We shall later come back to questions of 
ambiguities from the choice of the cut-off; for the time being, we shall adopt the same proper 
time regulator as in ref. p4|; then n((5^) has the form 



167r2 



dxx{l — x)T{0, xq) 



(3.3) 



where 



M2+Q2^(l-a;) 



and r(n, e) = / 



°° '^^e is the incomplete Gamma function. The 



X 

constituent quark mass Mq is the non-trivial solution {Mq ^ 0) to the gap equation 



Gs 



Mr 



exp - 




M. 



Q 



Ml 



A2 

X 



A2 

X J 



(3.4) 



The shape of the Adler function predicted by the ENJL Lagrangian is shown in Fig. 5 
(the dotted curve). As the figure clearly shows, ENJL does not interpolate correctly the 
intermediate region between the very long-distance regime and the short-distance regime; but 
what is more striking is that it also fails to describe hadronic data already at such low energies 
^■s ^ 4Mq {^Q — 300 MeV). The nai've expectation that this effective Lagrangian should 

The wanted pattern is the one which produces a sizable quark condensate. Other possible realizations of 
spontaneous chiral symmetry in QCD have been advocated by J. Stern and collaborators [^-|49[] and references 
therein. See however jl9|] . 
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be appropriate in describing the euclidean behaviour of Green's functions up to ~ A^, 
(since this is after all the scale that appears in the denominators of the composite operators 
in eq. (|3.1|),) has not been fullfilled |^. This failure can be delayed to slightly higher Q-values 
by introducing extra four-quark operators with higher derivatives in the lines proposed in 
refs. |]5T|, p^], but this introduces more unknown couplings and does not address the lack of 
confinement. 



0.04 




Q [ GeV ] 

Fig. 5 The Adler function predicted by the ENJL-Lagrangian in eq. J\ ) (dashed line) 
compared with the Adler function from the Aleph data (full line). 

What is the reason for this failure? In order to find guidance in answering this question, 
let us compare -4enjl(Q^) in sq. ( |3.2D with the expression which results from the LMD ap- 
proximation to QCD(oo) which we have been advocating in the previous sections, i.e. the 
approximation where 



-ImHc 

vr 



Nc 4. 



,{t) = 2f^M^5{t - M^) + —^-e{t - so) 



(3.5) 



167r2 3 

and where for simplicity we are switching off as corrections in the pQCD continuum term. 
This spectral function produces in the euclidean region the Adler function 

f2 T\/r2 \ AT A Q'i 



d 



+ 



N, 4 

167r2 3 So + 



(3.6) 



Clearly eq. ( |3.6| ) does not look like eq. (t3.2|). Even in the region of small we see that, 
firstly, eq. (|3.2D lacks Q'^/sq contributions that are present in eq. (^); and, secondly. 

This solves the "paradox" raised in ref. [Q. There is no overlap between the region where the ENJL- 
Lagrangian is operational and the region covered by the OPE. 

We would like to issue a warning message here: the continuum contributes to some of the 0{p^) (and 
higher) couplings of the chiral expansion. Consequently, saturating these coupling constants with only the 
resonance contribution, as is sometimes done by invoking e.g. VMD arguments, may be misleading. 
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also the dependence of the resonance contribution is misrepresented in eq. (O), at least 
insofar as n((5^) is dependent. All in all, the low-Q^ dependence of eq. (3]^) appears 
to be "incorrect" if expected to resemble the one in eq. ( |3.6| ) but, as we shall see, in a way 
that can be amended. The amendment will consist in adding the proper higher dimensional 
composite operators. We shall see that this addition is nothing but the implementation of 
the matching conditions inherent to any effective Lagrangian construction. 



3.2 Higher Dimension Operators and Confinement 

There will be two distinct classes of higher dimensional operators according to whether they 
originate from the continuum or from the prominent low energy narrow state, which we shall 
call for short the resonance. This is because chiral symmetry is realized differently in these 
two regions of the spectrum. In the case of the continuum, chiral symmetry is manifest by 
construction and the new higher-dimensional operators to be added to the Lagrangian in 
eq. ( |3.lD will contain solely combinations of external fields suppressed by the corresponding 
powers of sq, with chiral symmetry realized in a linear way. A possible example is the operator 

-tr{FL^,{x)aFj^''{x) + Fn^,{x)aF^''{x)} . (3.7) 
•So 

On the other hand, operators originating in the contribution from the lowest resonance must 
have spontaneous chiral symmetry breaking built in, and therefore they will appear as chirally 
symmetric in a nonlinear way; i.e. they will eventually generate effective couplings involving 
the Nambu-Goldstone bosons encoded in the [/-matrix and suppressed by a characteristic 
mass scale proportional to the QQ production threshold 2Mq. A possible example of an 
operator of this type is the dimension six operator 

^tr{Fz.^.(x)C/(x)t[aF]^^(:r)][/(x) +F^^,(:r)C/(x)[nFr(x)][/(x)t} . (3.8) 

In order to proceed to the explicit construction of these higher dimension operators it 
helps to have yet another perspective on the failure of eq. ( |3.1j ) to reproduce the expected 
structure of the LMD approximation to QCD(cxd) . When computing the imaginary part of 
IIenjl in eq. (3.2) one finds: 



1 , , ilmn(i) 

-Imn t ENJL = ; " , 2 , ; 72 • 3. 

(l-t^Ren(t))%(t^Imn(t)^^ 



The ENJL spectral function has a continuous shape, as shown in Fig. 6, whereas the large-iVc 
result should be a sum of narrow states; or at least a dominant narrow state as in eq. (|3.5| ). 
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0.6 0.65 0.7 

Vt [GeV] 

Fig. 6 The vector spectral function predicted by the ENJL-Lagrangian /|£ 

The shape in Fig. 6 seems reminiscent of the one corresponding to a "/?-resonance" with a 
"finite width" proportional to 

7 = t^^Imn(t) , (3.10) 

and, at first sight, it could even be taken as a success of ENJL. This shape however is more of 
a failure than a success because it involves a "finite width" which is proportional to Imll(t) 
and therefore of 0{N^) while in the large- A'^c limit of QCD, which after all >Cenjl in eq. ( |3.1D 
is supposed to approximate in some way, it should be of 0{1/Nc). The wrong large-A'^c 
behaviour of the ENJL width is actually due to the fact that the ENJL spectral function 
has on-shell QQ pairs (which involves an extra factor of Nc), i.e. precisely to the lack of 
confinement. What is needed here is to find the extra physics which generates the formal 
limit Imn(t) ^ in a consistent fashion, for then as 7 ^ 



7 



lim 

-° (i-t^Ren(t))- 



+ 7 



^Sil-t^^^Rem^ 



(3.11) 



and the narrow width limit would follow. 

The obvious way to perform this limit that comes to one's mind is to take the limit 
Mq — > CO. This is actually what happens in QCD(oo) in two dimensions 54 1, where the 
mass of the quark is pushed to infinity by the need to regulate the theory in the infrared 
region; and in two dimensions this is in fact how confinement takes place. Here this does not 
work because the constituent quark mass Mq given by the dynamics of /^enjl in eq. ( |3.lD 
is finite. One can nevertheless proceed as follows. Let us consider the very low-Q^ region. 
Since at = the operators of eq. (3J) produce the right physics, one can start at = 
and work the way up by examining what happens when trying an expansion in powers of 
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Q"^ /AMq. Because Mq is finite, the dependence of n(Q^) is a consequence of the fact 
that Imn(t) is nonvanishing. The dependence of tl{Q^) originates from the production 
of constituent QQ~y>^\ts in the physical spectral function; in other words from the lack of 
confinement. Indeed, the dependence is given by the dispersive integral 

n(Q2) = n(o) - q2 / imn(t), (3.12) 

J4M2 t(t + (y^j vr 

where ^lmn(i) is the spectral function associated with the constituent QQ^pairs. It is in fact 
this dispersion relation that holds the key to the remedy we are looking for. The crucial point 
is that eq. ( [j.l2D shows that t{[Q^) is analytic for < 4Mq. It can be expanded in powers of 
Q'^/AMq and therefore all the "wrong" dependence causing trouble can be systematically 



removed by adding to the Lagrangian in eq. {3A) an infinite set of local operators involving 
external fields of dimension n > 6: 

oo 

Aes = E CnOr- (3.13) 

n=6 

The lowest dimension is six because the first operator corresponds to the first derivative of 
the self-energy function n(Q^). The coefficients Cn in front of these local operators are known 
since they are adjusted so that all powers of Q'^/AMq coming from the Taylor expansion of 
n(Q^) cancel out in their contribution to the Adler function in eq. ( p. 21 ). For instance an 
explicit operator of dimension six which removes the "wrong" non-confining dependence of 

= n'(0) tv{uHD"Fj^'')U{DaFL^,,)} , (3.14) 
where as usual the trace is in flavour space, and n'(0) stands for 

The presence of r(l,e) in eq. ( p.lSj ) is due to the proper time regularization which we have 
adopted to define n(Q^) in eq. (|3.3D . The "wrong" non-confining dependence is regu- 
larization dependent and therefore, the counterterms which are added to cancel them will 
be regularization dependent as well. The net physical effect when the infinite set of local 
operators in eq. ( |3.13| ) is added, is formally the same as letting Mq — > oo in the ilmn(t) 
spectral function, namely n(Q^) n(0). 



3.3 Higher Dimension Operators and Matching Conditions 

The higher dimension operators in eq. ( [3.13 ) can also be viewed as part of the matching 



conditions which are required on the ENJL Lagrangian in eq. ( |3.1D to become an effective 
Lagrangian describing a narrow width confined vector state. Indeed, once these operators are 
included, effectively one is doing the following replacement in eq. 



nKNaL(Q') - n(o) ( -Q2__Vn(o)j , (3.16) 
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which is nothing but the expansion around = of 



and has now the form produced by an infinitely narrow state, i.e. 

dt 1 



noff(Q') = / — ^-Imn,a(t) , (3.18) 



with 

ilmn,«(t) = ^^S{t - A^cA^(87r2Gvn(0))-i) . (3.19) 
This is precisely the desired result. It leads to a hadronic contribution to the Adler function 

(iQ2 I M2 + Q2 



A.(Q^) = -Q'J^. f , (3.20) 



with 



VvM'^ = ^-T^ and = AT.A^ (8^2^yn(0))-^ (3.21) 



i.e., the same contribution as the LMD approximation to QCD(oo) formulated in eq. ( |3.6| ). 

At this point the reader may perhaps wonder about what all this improvement on ENJL 
has really brought us since it looks as if all we have accomplished is to reproduce eqs. ( |3.5| ) and 
( |3.6D , which after all we already knew from large- A'^c QCD arguments in the first place. The 
point, as we shall further develop in the next sections, is that imposing that the low-energy 
physics be governed by the effective Lagrangian in eq. ( p.l| ) with the infinite set of matching 
operators of higher dimension understood, constrains the parameters of the dominant low- 
energy states, which would otherwise be arbitrary. It relates, in particular, the lowest vector 
state to the lowest axial-vector state in a way which we will describe next. 



4 The Axial— Vector Channel 

We shall now be concerned with two-point functions 

KilU = i J d'^xe^'i-- < 0|r{A^(x)4(0)}|0 > (4.1) 
of axial-vector quark currents 

Al{x)=q{x)^^^^^q{x), (4.2) 

where A'^, as before in eq. ( |2.2[ ), are Gell-Mann matrices acting on the flavour triplet of u, d, s 
light quarks. In the chiral limit where the light quark masses are set to zero, these two-point 
functions depend only on one invariant function (Q^ = _g2 > q £qj, ^2 gpacelike) 

n^.(g)afe = (Qt^qu - 9f.uq'')IlA{Q^)Sab ■ (4.3) 
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The axial-vector two-point functions predicted by the ENJL Lagrangian have been discussed 
in ref. [24|. The self-energy function IIa{Q^) has the following form 

n^(Q2) ^ (4 4) 

where T[a{Q^) is now given by 

n^(Q2) = n(Q2) + Z!p!), (4.5) 

with n(Q2) the same function as the one defined in eq. (|3.3| ) and 

AT rl 

fiQ') = 8M^Q I dxTiO, xq) . (4.6) 

It is convenient to introduce the function 

^^^^ ^ ^ l + {Gv/AmM^J^dxr{0,XQ) ' ^^'^^ 



and then rewrite 11^ (Q ) as follows 



As in the case of the vector channel one notices that eq. ( |4.8| ) could match onto an axial Adler 
function 

A.(O')- ^ (2p 2flMl_\ AT, 4 

^^^^ ^--^dQ^y-Q^^ MfW ) + 1^3^^ ' ^^-^^ 

generated by a narrow width spectral function corresponding to the LMD approximation in 
QCD(oo) i.e., 

^lmn^(t) = 2f^6{t) + 2/lMl5(t - Ml) + r^p{t - ^o) , (4.10) 

provided that both tl{Q^) and f'^{Q'^) get their dependence, induced by QQ discontinu- 
ities, canceled by the inclusion of appropriate higher dimensional composite operators in the 
external fields. 

We emphasize that the procedure of introducing an infinite number of local operators 
adds no unknown parameters, as the coefficients accompanying the operators are fixed by the 
coefficients of the Taylor expansion at zero momentum of the functions Il{Q^) and /"^{Q"^) 
which are explicitly known. Furthermore, one also has to add the operators that fix the 
contribution of the continuum, as was explained in the case of the vector channel in subsec. 3.2. 
Once this is implemented the masses and coupling constants of the dominant low-energy 
states in the vector and axial-vector channels are completely fixed in terms of the parameters 
Gs, Gv and A,^ in a way which we summarize below. 
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i) The gap equation. 

The gap equation in eq. ( |3.4| ) relates Gs to the ratio of the constituent quark mass Mq 
and the scale 

1 

— = er(-l,e), with (4.11) 

ii) The coupling constant qa- 

The coupling constant gA is governed by the size of the Gy coupling 

gA = gA{0) = (1 + 4Gyer(0, e))-^ . (4.12) 

iii) Masses. 

The masses of the vector and axial-vector states are then fixed as follows 

Mv = I^T^,=^K-r^^ Ml = ^. (4.13) 

2 Gv r(0, e) ^l- gA gA 

iv) The Coupling Constants. 

The coupling constants of the lowest narrow states are then also fixed: 

/' = ^4M25^r(0,e), (4.14) 
N 2 

/I = Si/f. (4.16) 

With these results, the 1st and 2nd Weinberg sum rules, which in the case of LMD imply the 
well known relations 

f^ + flMl = flMl, (4.17) 
flM\ = f^M^, (4.18) 

are then automatically satisfied. It appears then, that just removing the non-confining terms 
in the vector and axial-vector two-point functions predicted by the usual ENJL-Lagrangian 
guarantees their correct matching to the QCD short-distance behaviour. This shows, in ret- 
rospect, that the choice of an initial ENJL-ansate as a starting effective Lagrangian was quite 
a good one. However, as we shall later see, the removal of the non-confining contributions 
in Green's functions beyond two-point functions is not enough to guarantee, in general, the 
correct matching to the QCD short-distance behaviour and further local operators have to 
be introduced. 



4.1 The Electromagnetic 7r+ — n'^ Mass-Difference 

In the chiral limit, and in the presence of electromagnetic interactions to lowest order in the 
fine structure constant a, the hadronic spectral functions in eqs. (|]5D and ( CT ) induce an 
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electromagnetic mass for the charged pfon, which is finite when the spectral functions obey 
the 1st and 2nd Weinberg sum rules ^ This results in the expression 



a 3 



,r+ |em 



Ml 



log 



Ml 

M2 



which in terms of the parameters Mq and qa defined above reads 



M^ — 

7r2 '3(1 



QA 



9 A) 



■log — . 
gA 



(4.19) 



(4.20) 



We show below in Fig. 7 a two-dimensional plot of the resulting electromagnetic pion mass 
difference Am^ = — m^o as a function of Mq and when these parameters are varied 
in the ranges 200 MeV < Mq < 400 MeV and 0.4 < gA < 0.6. Experimentally [||], Am^r = 
(4.5936 lb 0.0005) MeV. This corresponds to the plateau shown in Fig. 7. The intersection of 
eq. ( [4.20 ) with this plateau results in a line which correlates Mq and gA- The figure shows that 
a very reasonable choice of values for these two parameters can reproduce the experimental 
result. We shall be more precise about the values of these parameters in subsection 5.2 



Am 




350' 



Fig. 7 The electromagnetic pion mass difference Atti-tt in MeV versus gA and Mq (MeV) 



5 The Low— Energy Effective Lagrangian 

The degrees of freedom of the resulting effective low-energy Lagrangian are a pseudoscalar 
Goldstone field matrix, a vector field matrix V{x), a scalar field matrix S{x), and an axial- 
vector field matrix A(x). They correspond, therefore, to the wanted LMD spectrum of 0^, 1~, 

For an introduction to this subject, see e.g. the lectures in ref. where references to the earher work 
are also given. 
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0"*", and 1"*" hadronic states. We postpone the technical discussion concerning the removal of 
QQ discontinuities in Green's functions with external scalar sources to another paper. Here, 
we shall limit ourselves to quote the results we obtain for the effective 0{p'^) chiral Lagrangian 
of the Goldstone modes after integrating out the physical scalar fields. It is this dynamical 
effect which produces the L5 and Lg couplings as well as part of the L3 coupling which we 
give below. The remaining interactions between the pseudoscalar Goldstone field matrix and 
the vector and axial-vector field matrices produce the rest of the couplings in the following 
effective Lagrangian, valid up to 0{p'^) and also for some of the higher order terms, 



--efl 



J/2 [tr (L»^C/D^C/t) + tr (x?7^ + U^x)] + Ctr(Qj?C/QiC/t) 

tr [V^uV^"" - 2M^Vf,V>'] 



4 

-itr \V....V^"' - 9/\/2 

-Jtr[V^'^"-2Ml^^^''] 
1 



[/ytr ( V/f ) + i~gvtT {Vf,,[e, el) + Atr (V/-")] 



2V2 

+Li (tiD^U^D^'Uy + 2Litr (^D^U^ D^Utv Di'U^ D^u) 
+L3tr {p^U^D^'UDM^D''u) 

+L5tr [L>^C/tDMC/ {x^U + C/tx)] + Lgtr [x^Ux^U + xU'^xU'^) 

-iZgtr (f^^D^UD.U^ + F^^D^U^D^U) + Liotr {u^ F^'^UFl^.u) 

+HiiT {F^'^Fr^, + F^Fl/..) + H2tv (x^x) ■ (5.1) 



Here, 



yij-i' ^ y/xyv _ v^F'', Ai^" = - VA^" ; (5.2) 

with V the covariant derivative defined as 



v^i? = a^E + [r^,i?], (5.3) 



and the connection 



We recall that and are external vector and axial-vector matrix field sources, (not to be 
confused with the hadronic matrix fields and Af^,) and Fj^'^ and Fj^'^ are the external field 
strength matrix tensors associated with the Ifj, = — and = + sources. We are 
also using the notation U = with 

/r = Ci^rC^±e^i^re and C^ = ^^^D,U^K (5.5) 

This is the effective Lagrangian which results from the ENJL-ansate after removal of non- 
confining QQ discontinuities. As already mentioned, this guarantees the correct matching to 
the QCD short-distance behaviour for two-point functions only. The couplings gv, Lg, and 
Li , L3 are correlated with three- and four-point functions which still require the introduction 
of extra higher dimension operators to guarantee the correct matching to the leading QCD 
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short-distance behaviour. This is why we put a tilde on them, and we shall discuss how they 
have to be modified in the next subsection. All the constants appearing in the Lagrangian 
in eq. ( |5.1| ) are functions of Gs, Gy, and alone or equivalently g^, Mq and Tq = T{0, e). 
The masses My and have already been given in eq. ( 4.13| ), and the couplings fy and 
in eqs. ( |4.14| ), (4.15) and (4.16). The constant G is related to the electromagnetic vr"*" — vr'^ 
mass difference in eq. ( [4.20| ), 



2e^G 

J IT 



(5.6) 



One could in principle consider as well terms of O(e^p^) but this goes beyond the scope of 
this paper. At this stage i.e., after removing QQ discontinuities in the initial ENJL-ansate, 
the result for the VPP coupling gy is 



~9l 



Nr 1 



1 



167r2 6 

and the results for the other couplings in eq.(|5.lD are: 



Li 



is 



Nc 1 
167r2 48 
Nr 1 r 



167r2 



(i-5i)'ro, 

(1 - al? + 2g\ 



Nr 1 



L 



10 



167r2 4^ 

16^16^^^° 

N^l 
167r2 6 ^ 

Nc -1 
167r2 6 



(i-9i)ro 



(5.7) 

(5.8) 
(5.9) 
(5.10) 
(5.11) 
(5.12) 
(5.13) 



Integrating out the vector and axial- vector fields in the Lagrangian of eq. (5.1) produces 
terms of 0{p^). However, to that order, there are two extra sources of other possible con- 
tributions which must be taken into account. One is the effect of integrating out the vector 
and axial- vector fields in the other 0{p^) interaction terms which are not included in C^g 



in eq. (|5.1|) , (terms like e.g. e/ii/p(Ttr(y^,^'^,^'',^'^);) the other source, as already mentioned, is 
the effect of new operators coming from the matching of Green's functions beyond two-point 
functions to short-distance QCD. 

The Adler two-point functions which we have discussed in the previous sections, i.e. the 
expressions in eqs. (|3.6D and (|4.9| ) with the masses and couplings as given in eqs. ( 4.13|) and 



( 4.14 ), ( 4.15 ), ( 4.16 ), are the result of the full resummation of all the relevant terms in the 
effective Lagrangian. Expanding these Adler functions in powers of is, of course, another 
way of obtaining the couplings of the corresponding operators in the effective Lagrangian 
which govern the behaviour of these functions. 



5.1 Dynamical Symmetries of the Low— Energy Effective Lagrangian 

The results we find at this level for the low-energy couplings show already a lot of interesting 
symmetries. It is worthwhile describing them and comparing them with the dynamical rela- 
tions discussed in ref. |27|. For the sake of comparison it is more convenient to rewrite the 
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results we have found for the couphngs /y, /a and gy in the following way 



fv 

fA 

9v 



1 



1 



9 A] 



9 A My ' 



1 



1+9A 



9A- 



2 " My 

The results for the couplings in eqs. (5.8) to (5.13) can then be written as follows: 

= 1 



o9v 



-6L1 + qaL^ , 



Nc 1 
167r2 16 

^fvgv , 



9a^o 



fl) 



(5.14) 

(5.15) 
(5.16) 

(5.17) 

(5.18) 
(5.19) 

(5.20) 
(5.21) 
(5.22) 



Some of the relations encoded in the Lagrangian (|5.1|) turn out to be common to those 



discussed in ref. |27] in connection with the formulation of chiral Lagrangians of vector and 
axial-vector mesons obeying short-distance QCD constraints. They are the following: 



The 1st and 2nd Weinberg sum rules in eqs. (|4T^ ) and (|I1). 

The relationship between the constant Liq and the /y and /a couplings in eq.( |5.22| ). 

The relationship between the constant Zg and the /y and gy couplings in eq. (|5.2l| ). 

The relations: ^ 

Li = —gy , and Lalvector = — 6L1 , (5.23) 



where Z/3I vector means the contribution induced by vector exchange to L3. In ref. [^] it is 
shown that these two relations follow from the assumption that the forward scattering 



amplitudes of two pseudoscalar mesons obey the Froissart bound |56 



The fact that the effective Lagrangian in ( |5.1| ) obeys automatically these relations shows again 
that the ENJL-Lagrangian ansatz plus the corresponding set of local operators required 
to remove the non-confining QQ-discontinuities which it generates is already a very good 
effective Lagrangian to describe the LMD approximation to QCD(cxd) . There are however 



other relations in ref. |27] which the effective Lagrangian in (|5.1| ) does not quite satisfy. Let 
us describe them. 

i) In ref. it is shown that the assumption that the pion electromagnetic form factor F{(p') 
with the charge normalization condition F{Q) = 1, obeys an unsubtracted dispersion 
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relation, and the saturation of the dispersive integral with just the p-pole leads to the 
relation 

fv9v^ = l, (5.24) 
with gy the physical /ovrvr coupling constant. 

ii) A similar assumption of an unsubtracted dispersion relation for the axial form factor in 
vr — > efe7, with the same saturation assumption of the dispersive integral by a resonance 
pole was shown in |27] to lead to the relation 



fv = 2gv , (5.25) 



which is sometimes called the KSFR-relation [58,5 
The combination of i) and ii) results in 



fv = 2gv = . (5.26) 



The combination of i), ii) and the two Weinberg sum rules in eqs. (P^) and (p^ imply 



Ma = V2Mv (5.27) 

and therefore 



'i) =-5Mjlog2. (5.28) 

/ EM vr Z 



With and My as independent variables, the other couplings are then all fixed: 



When compared to the results in eqs. (|5ll , (|1|) and (|5l5|) it looks as if there IS no unique 
solution for gA which can reproduce all these constraints simultaneously. The value gA = \ 
reproduces the result Ma = \piMy and the couplings for /y and jA in eq. ( |5.2g|) , but not gy . 
On the other hand, the electromagnetic pion form factor constraint in eq. ( |5.24D can only be 
satisfied '\i gA = Ij while the axial-form factor constraint in eq. ( p.25| ) requires = 0; but 
both these two extreme values are excluded by general large- A'c arguments [^]. 

From the description above, there follows that the question of whether or not the electro- 
magnetic pion form factor and the axial form factor in vr — > ez^e7 obey unsubtracted dispersion 
relations is an important issue. In QCD it is indeed expected that these form factors do obey 
unsubtracted dispersion relations, and there is no reason that the restriction of full QCD to 
QCD(oo) with LMD spoils this property Let us then consider the electromagnetic pion 
form factor F(C^) in somewhat more detail. The effective Lagrangian in eq. (|5.l| ) gives: 



Short-distance QCD arguments for a, fall-off of the electromagnetic pion form factor can be found 

in ref. The fall-off also follows rigorously from an analysis of QCD short-distance constraints 

on three-point functions saturated by a vector state and pseudoscalar Goldstones only. 
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The first term 1 is the 0{p^) contribution; the second term comes from the Lg-coupling; the 
last one from integrating out the vector field V from the 0{p^) interaction terms. Recall that 
Lg = ^fvgv, with fv and gy defined in eqs. (5.13) and (5.14). Clearly, 

M2 Q2 / q2 



i + gA 



(5.32) 



The problem with this form factor is that it does not satisfy the QCD requirement that 
asymptotically in it should fall as It does not have either the expected LMD behaviour 

Notice that here we are talking about an on-shell form factor (the pions are on-shell), and its 
imaginary part is supposed to have one state only, the lowest vector meson state. There is no 
perturbative on-shell continuum here. In order to implement the correct large-Q^ behaviour 
we are then forced -once more- to introduce a set of an infinite number of local operators in 
external fields. These are operators of the second class type discussed in section 3.2 . Their 
contribution to F{Q'^), let us call it f{Q'^), when re-summed must give: 



2 

for then, the overall form factor will be 



Fm - 1 - ^+3A _ i-gA _ 



as wanted. The overall effect of this infinite number of local operators is a simultaneous 
modification of the coupling constant gy and the constant Lg in the effective Lagrangian in 
(|5.l|) as follows: 

1 + ffA n fn fl + gAl-gA\ n U ,^ „„^ 

= VI -5^^ = gv , (5.37) 

Lg = -^-(1 - gDTo -^-(1 - gA)To = Lg . (5.38) 
^ levr^e^ 167r2 3^ wa; u y \ j 



and 



A similar modification of the coupling constants Li, L3 also has to be made, accordingly, so as 
to preserve the relations in eq. ( 5.23| ) which ensure the short-distance matching for four-point 



functions, with the result 



2N2p 1 
16^2 48 167r2l2' 

L3 = -6L1 + gAL5 -6L1 + 5AL5 = L3 . (5.40) 



Li = T7^j^ii- g\fro ^T^-^ii- gAfTo ^ l, , (5.39) 
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The new coupling constants preserve the symmetries 



L3 -- 

Lg -- 

and, now, they also satisfy the relation 



L 2 



1 



-6L1 + L5 
-J'vgv , 



fvgv 



M2 



(5.41) 
(5.42) 
(5.43) 



(5.44) 



Let us next consider the axial form factor Gyi(Q^) in vr — > ei"y. Again, the condition 
of saturation of the imaginary part with one axial-vector state, and no subtractions in the 
dispersion relation implies |27]: 2fvgv — /y = • It turns out that now, using eqs. ( 5.14 ) and 
(5.44), there is indeed a solution which satisfies this constraint, namely: 



(5.45) 



The net effect is that the number of independent input parameters has now been reduced 
from three to two. 



5.2 The LMD Effective Lagrangian of Large-iV^ QCD 



We finally find that after correcting the initial ENJL-Lagrangian ansatz for the unwanted 
non-confining QQ-discontinuities and demanding the correct matching between the leading 
QCD short-distance behaviour of the VPP and PVA three-point functions and the asymp- 
totic behaviour of the corresponding hadronic form factors, the resulting low-energy effective 
Lagrangian which describes LMD in QCD(oo) , to 0{p^) in the chiral expansion, has the 
simple form 

= [tr (^D^UDf'U^) + tr {xU^ + U^x)] + Cti{QRUQLU^) 



tr [V^.^V^'' - 2M^V^V^] 



--tr [A^.A'^'' - 4M^A^A^'] 



1 fn 



4 My 



tr (2V/f + iV^u[e,e] + A^^f: 



+Li ftr DaU^D^'U] + 2Litr (dM^ D^UiiD^'U^ D'^U 



+L3tr (dM^D^'UD^U^D'^U 



+L5tr \D^U^D^'U ( x^U + U^x)] + ^str f x^Ux^U + xU^xU^ 



-iLgtr [f^'D^UDM^ + F^'D^U^D.U) + Liotr (u^F^'UFl 
-Hit! {F^'Fr^, + F^'Fl^,.) + i72tr (x^x 



(5.46) 
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with, so far, two free parameters Mq and Tq which are such that: 



f^ = J^2M^To and M^ = ^mI = 6M^. (5.47) 



167r2 

The constant C is then fixed, 

^ = ^9M^log2, (5.48) 
and the Gasser-Leutwyler Lj constants are ah proportional to -^^Tq, i.e. to 3-^2-: 

Q O TV 1 

6Li = 3L2 = — = 4L5 = SLg = -Lg = -^lO = TTT^ • (5.49) 

Notice that the dependence on the regularization chosen in the ENJL-ansate only appears as 
an undetermined overall constant Tq. Different regularizations will give different "expressions" 
for Tq, but once the value of Tq is traded for a physical observable, all the predictions are 
regularization independent. 

When restricted to the vector and axial-vector sector, the Lagrangian in ( ^.461 ) is entirely 



analogous to the phenomenological Lagrangians discussed in ref. |27]. They are just different 
formulations of the same physics which describes the LMD of the 0~, 1~ and 1^ hadronic 
states in QCD(oo) . The predictive power of the low-energy effective Lagrangian in 



as compared to various phenomenologically inspired Lagrangian formulations of resonances 



discussed e.g. in refs. [pq-31| and references therein lies in the fact that in our case the vector, 
axial-vector and scalar couplings are all correlated by one mass scale Mq and one dimen- 
sionless parameter Tq, which can be traded for /^r and My- In order to predict the 0{p^) 
Lj-constants e.g., the phenomenological approaches require f-,^. My, the scalar couplings 
Cm, Cd, and Ms as a minimum of input parameters. Going from 0{p'^) to 0{p^) in a phe- 
nomenological approach requires many more new input parameters, while in our case they 
are all fixed by the same Mq and Tq parameters, provided that the matching constraints to 
short-distance QCD for the new Green's functions which will appear at 0{p^) can be made 
consistently with the LMD 0~ , 1~ , O"*" , 1"'' spectrum only. In fact, we have seen that at O(p^), 
the matching constraints to short-distance QCD beyond two-point functions have actually 
reduced the number of three free parameters in the initial ENJL-Lagrangian ansatz to two. 
We have not yet made an exhaustive study of this important question beyond O^p"^). 

One way to quantify the predictive power of the LMD approximation to QCD(oo) which 

we have been considering is to fix the remaining two parameters Mq and e = ( or Tq) 

from a least square fit to experimentally well known quantities and compare them with the 
resulting predictions. The choice of low-energy constants is the one in the first column of 
Table 1. 
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Table 1. Low-Energy Constants, Experimental Values and Least Square Fit. 



Parameter 


Equation 


Exp. Value 


Fit Error 


Predicted Value 


T 


(p.49 


) 


(1.4 ± 0.3) X 10 


±0.3 X 10 


1.7 X 10~3 


L3 


(M9 




(-3.5 ± 1.1) X 10-3 


±1.1 X 10-3 


-4.4 X 10-3 


L5 


(K49 




(1.4 ±0.5) X 10^3 


±0.5 X 10-3 


1.3 X 10-3 


Ls 


(M9 




(0.9 ±0.3) X 10^3 


±0.3 X 10-3 


0.64 X 10-3 


L9 


(^.49 




(6.78 ±0.15) X 10-3 


±0.15 X 10-3 


6.8 X 10-3 


Lio 


(M9 




(-5.13 ±0.19) X 10-3 


±0.19 X 10-3 


-5.1 X 10-3 




(^.47 


) 


(92.4 ± 0.3) MeV 


±9 MeV 


(87 ±3.5) MeV 




(^.48 


) 


(4.5936 ± 0.0005) MeV 


±0.5 MeV 


(4.9 ± 0.4) MeV 


Mv 


(B.47 




(768.5 ±0.6) MeV 


±75 MeV 


(748 ± 29) MeV 


Ma 


(M7 




(1230 ±40) MeV 


±200 MeV 


(1058 ± 42) MeV 



The equations which define them in the text are given in the second column. The third 
column quotes their present experimental values with their errors ^ The fourth column 
shows the errors which we have taken for the fit. The errors for the input parameters with 
a mass dimension are taken larger than the corresponding experimental error, to allow for 
possible effects due to next-to-leading l/A^c~corrections and chiral corrections. From the nine 

constraints in Table 1 and with the two parameters Mq and e = as free parameters we 

obtain a fit with a = 4.2. The corresponding values of the fit parameters are then 

Mq = (305 ± 12) MeV , and e = 0.0706 ± 0.0031 . (5.50) 

The resulting values for the Lj-constants and the other parameters are given in the last 
column. Typically, the errors for the predicted Lj's are less than 2%; but of course this 
ignores the systematic errors of the approximation which we are adopting. Figure 8 below 
shows a plot of the predicted L[s in eq. ( |5.49|) versus Fq. 

^^We use the recent determination of Lio and the value of Lg quoted in ref. p%- 
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0.5 1 1.5 2 2.5 

r(0,8) 



Fig. 8 Plot of the predicted Li constants in eq. ^5.4^ versus Tq. The points with error 
bars are the present phenomenological determination of these constants from experiment. 

It is remarkable how well the phenomenological determination of these constants intercept 
the predicted linear dependence at the same value of Fq. This is a clear improvement, both 
in the quality of the fit (with one less parameter!) and in conceptual simplicity, to previous 
fits made with the ENJL-Lagrangian [^,^ which were already quite good. 

All the other parameters which have appeared at various stages of the discussion in the 
text are now fixed. For example the couplings of the initial ENJL-Lagrangian are found to 



28 



be 

G5 = 1.28 ±0.01, = 1.65 ±0.04, and = (1150 ± 52) MeV . (5.51) 

Another parameter which has been discussed in the text is the onset of the QCD continuum 



So in the vector and axial-vector spectral functions in eqs. ( |3.5[ ) and (4.10). This parameter 
was fixed by the requirement that in the OPE of the corresponding two-point functions there 
is no local operator of dimension d = 2 in the chiral limit. This leads to eq. ( 2.21[ ) and hence, 



using the determination for fyMy which follows from the effective Lagrangian in ( |5.46| ) fixes 
the value of sq to 

so(l + ...) = 6M^ro~ 1.20GeV^ (5.52) 

where for simplicity we have neglected gluonic perturbative corrections. Once the value of sq 
is fixed, the matching to higher 1/Q^-powers in the OPE fixes the corresponding combinations 
of condensates of gauge invariant local operators. For example, matching the ^ powers in 
the Adler function in the chiral limit, constrains the gluon condensate to satisfy the sum rule 

— {G^'^Gl^) = -Af^M^ + ^\sl{l + ■■■). (5.53) 

which in terms of the two parameters Mq and Fq and with neglect of pQCD corrections 
results in 

f (GrC^J ^ ^{2V^Mq)^To (Fo - 2) . (5.54) 
Numerically, this corresponds to a value 

«s(GrG'^.>=^ 0.048 GeV^ (5.55) 

which is quite compatible with the phenomenological determinations (see eq. ( |2.14| )). 

Finally, the resulting shape of the Adler function, which was the starting point of our 
analysis after all, is the dashed curve shown in Fig. 9 below. In this curve, the gluonic cor- 
rections in the matching with the pQCD continuum have been included. The agreement with 
the experimental shape (the full line) is not as good as the VMD one with the experimental 
values of Mp and fp plotted in Fig. 4, but it is not bad. 
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Fig. 9 The resulting Adler function of the LMD approximation to large-N^ QCD (dashed 
curve) compared to the experimental data (full line). 

As discussed in ref. [19|, in the LMD approximation to large- A'c QCD, and for given values 
of the masses of the lowest vector and axial-vector states, all the local order parameters which 
appear in the OPE of the ^lr{Q^) correlation function are fixed by the meson masses; in 
particular the matching of the powers leads to the sum rule 



vr 



which from eqs. ( |5.47 ) can also be written as follows 



(5.56) 



(5.57) 



Using the values for the parameters Mq and e which result from the least square fit, i.e. the 
values in eq. ( p. 50 ), and the value of at the matching scale in eq. (5.52) we obtain 



{iP^) = [(-288 ± 12) MeV]' 



(5.58) 



on the high ball park of the phenomenological determinations (see eq. ( |2.15| )) but certainly 
not an unreasonable value |^. 

One of the reasons why the ENJL-Lagrangian ansatz was chosen in the first place is that 
above the critical value Gs > 1, it produces a quark condensate, which in terms of Mq and 
e is given by |^] 



167r2 



4MSr(-l,e). 



(5.59) 



^®In fact, the scale for th e qua rk-condensate in eq. (5.59) is so — 1.2 GeV^, which is larger than the 1 GeV'^ 
reference scale used in eq. (2.15). The comparison at the same SQ-scale is slightly better. 
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So far, however, we have not used this equation at ah. // one is wilhng to identify this 
quark condensate to the one of QCD(c!o) at a scale /i ~ A^; an interesting relation between 
the pQCD coupling constant and non-perturbative parameters appears then when inserting 



eq. ( 5.59 ) in the r.h.s. of eq. ( 5.57 ) and using eq. ( 4.11 ) as well: 



asNc 



9G| 



TT 



4 (7? / 



(5.60) 



This relation can be viewed as a matching constraint between pQCD above the scale A-,^- and 
the non-perturbative effective realization in terms of hadronic degrees of freedom below A^. 

[Notice that T (o,M^/A'^^ = log{A'^/M^) - 7^ H .] Numerically, the central value of the 

r.h.s. of this equation which follows from the fit values in eqs. ( |5.50 ) and (5.51) is 0.63, while 

0.48, showing that the matching, at least in the vector and axial- 



gsNc 



lllog(A^/A: 



MS 



M2 



vector channels, is not unreasonable. In fact, one could use eq. ( ^.60|) to determine e 
in which case the low-energy parameters would all be constrained by just one mass scale. 



6 Conclusions and Outlook 

We have shown that there exists a useful effective Lagrangian description of a well defined 
approximation to QCD in the large-A'^c limit. This is the approximation which restricts the 
hadronic spectrum in the channels with quantum numbers 0~, 1~, 0^ and 1+ to the lowest 
energy state and treats the rest of the narrow states as a QCD(oo) perturbative continuum, 
the onset of the continuum being fixed by consistency constraints from the operator product 
expansion. The degrees of freedom in the effective Lagrangian are then a nonet of pseudoscalar 
Goldstone particles which are collected in a unitary matrix U{x), and nonets of vector fields 
V{x), scalar fields S{x) and axial-vector fields A{x) associated with the lowest energy states 
of the hadronic spectrum which are retained. We have derived the effective Lagrangian by 
implementing successive requirements on an ENJL-type Lagrangian (see eq. (^3)) which we 
have chosen as the initial ansatz. The first requirement is to eliminate the effects of non- 
confining QQ discontinuities by introducing an infinite number of appropriate local operators 
with couplings which can be fixed in terms of the three parameters of the starting ENJL- 
Lagrangian itself, i.e. Gs, Gy and the scale A^. We have shown that the matching of 
the two-point functions of this effective Lagrangian to their QCD short-distance behaviour 
can be systematically implemented. In particular, the 1st and 2nd Weinberg sum rules are 
automatically satisfied. The resulting Adler function shown in Fig. 9, when confronted with 
the experimentally known curve, shows a rather good interpolation of the intermediate region 
between the two asymptotic regimes described by perturbative QCD (for the very short- 
distances) and by chiral perturbation theory (for the very long-distances). 

For Green's functions beyond two-point functions, the removal of the non-confining QQ 
discontinuities produced by the initial ENJL ansatz is however not enough to guarantee in 
general the correct matching to the leading QCD short-distance behaviour and further local 
operators have to be included. We have discussed this explicitly in the case of the VPP 
and VPA three-point functions, and shown that the matching with the QCD short-distance 
leading behaviour which follows from the OPE restricts the initial three free parameters of the 
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ENJL-Lagrangian ansatz to just one mass scale Mq and a dimensionless constant Mg/A^. 
The resulting low-energy Lagrangian in the vector and axial-vector sector, and to 0{jp'^) in 
the chiral expansion, coincides with the class of phenomenological Lagrangians discussed in 
ref. which also have two free parameters and /^/M^. In this respect, this explains the 
relation to QCD of the phenomenological VMD Lagrangians discussed in ref. |27|: to 0{p^), 
they can be viewed as the effective low-energy Lagrangians of the LMD approximation to 
QCD(oo) . On the other hand, the fact that the resulting low-energy Lagrangian coincides 
with the phenomenological VMD Lagrangians discussed in ref. demystifies to a large 
extent the role of the ENJL-Lagrangian itself as a fundamental step in deriving the low- 
energy effective Lagrangian of QCD. The ENJL-Lagrangian turns out to be already a very 
good ansatz to describe in terms of quark fields degrees of freedom the LMD approximation 
to QCD(oo) and this is why it is already quite successful at the phenomenological level; but, 
as we have seen, when the non-confining QQ discontinuities are systematically removed, the 
phenomenological predictions improve even more. 

There is an advantage, however, in starting with the ENJL-Lagrangian as an ansatz, and 
that is that in this description of the LMD approximation to QCD(oo) , all the couplings to all 
orders in xPT are clearly correlated to the same two free parameters. For example, as shown in 
eqs.( 5.4E ), the -L5 and Lg constants, as well as part of the contribution to the L3 constant which 
result from scalar exchanges, are now proportional to the same parameter Fq, while in a purely 
phenomenological description in terms of chiral effective Lagrangians which include resonances 
as discussed e.g. in refs. p8|-^ and references therein, these low-energy constants require new 
phenomenological input. It is a very interesting question to ask whether or not a systematic 
implementation of short-distance OPE constraints on phenomenological Lagrangians leads 
to the same description, to all orders in the chiral expansion, as the one starting with the 
ENJL-Lagrangian ansatz. As we have said already, we have not yet undertaken this study in 
an exhaustive way. 

Finally, we wish to insist on the phenomenological successes of the LMD approximation 
to QCD(cxd) , as demonstrated by the results collected in Table 1 and the plots in Figs. 8 
and 9. These results show that this is indeed a very good approximation to full fledged 
QCD. It seems now worthwhile to apply it to the calculation of couplings of 0{p^) and 
0{e'^p'^) as well. One can also reconsider non-leptonic weak interactions in the light of this 
effective Lagrangian framework, where we expect a much better matching between the long- 
distance evaluation of matrix elements of four-quark operators and the short-distance pQCD 
logarithmic dependence of the Wilson coefficients than e.g. in the cut-off calculations of 
refs. p^,|6^] or the ENJL-calculations reported in ref. 
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